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INTRODUCTION 

In January of 1975 work was started at Georgia Tech on a project 
aimed at gaining a better understanding of the various design parameters 
which, affect steady state and transient operation of cryogenic heat 
pipes. This report briefly describes the progress made on the project 
during the period January through December of 1975. Financial support 
has come from NASA under grant NSG-2054 and the work has been monitored 
by Jack Kirkpatrick of Ames Research Center and Stan Ollendorf of 
Goddard Space Flight Center. One M.S. thesis (reference 1) which is 
directly related to the project was published in June of 1975 and a 
second thesis in the area is currently being prepared and should be 
published about July 1976. It is anticipated that several papers will 
be published in recognized technical journals over the next few years 
as a result of the work. 

Heat Pipe Under Study 

A 304 stainless steel heat pipe with slab type capillary structure 
and nitrogen as the working fluid was studied in the temperature range 
of 60^K to 120“K. The pipe is 1.27 cm in outside diameter and 9.14 cm 
in total length. Figures 1 and 2 show geometry of the pipe and the 
configuration of the capillary structure. In the transient studies, 
described in detail in this report, saddles are included at evaporator 
and condenser ends and a radiator is included at the condenser end. 
Summary of Results to Date 

The work performed under the grant is divided into two main areas. 
The first area includes development of accurate steady state equations 
for predicting capillary limitations and development of equations for 




End Cap (2) 

Figure 1- General Layout o£ Heat Pipe 
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predicting thermal resistances. Extensive computer surveys were run 
using the equations developed and it was found that a heat pipe of the 
type described above could transfer as much as 30 to 40 watts with an 
overall temperature drop of a few degrees Kelvin. The second area of 
study was the transient operation of cryogenic heat pipes. The problem 
was first studied with the aid of an analog computer and currently a 
digital computer approach is being pursued. Results of the analog 
work, which is limited to small transients and does not account for 
fluid dynamics, show the effects of small step temperature changes and 
the effects of fast and slow sine wave variations of small amplitude 
on heat pipe operating parameters. The digital program now under 
development will be much more powerful in that large transients, even 
including the case of start-up from the supercritical state, can be 
studied. At present a simplified model, which does not include fluid 
dynamic effects, is being successfully examined on a digital computer. 

The next step is to incorporate fluid dynamic effects. 

Future Plans 

The main theoretical effort is now being directed towards making 
the transient model, which is being examined on a digital computer, more 
realistic and more flexible. It is expected that this work will require 
another six months assuming that current levels of effort are maintained. 

Some effort is now being directed to planning some low temperature 
heat pipe experiments which could be used to generate data for verifying 
and modifying both steady state and transient models. In the near term 
the experiments would be carried out in a laboratory on earth. Also 
some initial thought has been given to the idea of a space experiment. 
The National Aeronautics and Space Administration currently is designing 
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a "Long Duration Exposure Facility" and a "Sky Lab" either of which 
would be ideal for this type of space experiment. 
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RESULTS TO DATE 

During Calendar 1975 significant progress has been made both in 
the steady state and transient parts of the work. Computer programs 
are now on hand for predicting capillary limitations, and steady state 
thermal resistances for slab type cryogenic heat pipes. An analog 
scheme has been developed which handles small transients and work is 
progressing on a much more powerful digital solution scheme which will 
be able to handle large transients including startup from the super- 
critical region. Each of these steady state and transient approaches 
will be discussed in detail. 

Steady State Thermal Resistances 

In the typical design of a heat pipe, little attention is given to 
predicting thermal resistances. This is the case because accurate pre- 
dictions are extremely difficult in most systems. It is not .uncommon 
to underestimate overall heat pipe temperature drops by an order of 
magnitude. 

In the present analysis resistances are considered in the pipe 
wall at the condenser and evaporator, in the layers of capillary material 
around the circumference of the evaporator and condenser surfaces, in 
the fluid gaps between layers of capillary material, at liquid vapor 
interfaces, and in the vapor region. In addition it has been asstmed 
that the circumferential portion of the capillary structure partially 
drys as heat transfer is Increased towards the capillary limitation. 

This problem is discussed in detail in references 2 through 5. Results 
of several studies, see reference 2, suggest that in long heat pipes' 
part of the condenser surface may not he active at relatively low heat 
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transfer rates. This effect has been considered in developing a thermal 
resistance model at the condenser end. 

The following nomenclature is used in computing thermal resistances, 
d^ distance between centers of screen filaments 
g conversion factor 

enthalpy of vaporization of 

conductivity of stainless steel screen filaments or of the 
stainless steel pipe 

k» conductivity of liquid nitrogen 

k^ conductivity of stainless steel pipe 

k^ conductivity of the liquid filled screen portion of the wick 

% axial length 

ax 

£ active condenser length 

C^3L 

length of condenser at design conditions 

£ length of evaporator section 

9 . effective length of vapor path 

efi 

N number of screen layers 

p^ vapor pressure 

Q heat transfer rate 

^MAX limitation 

R ideal gas constant 

r^ outer radius of pipe wall 

r^ inner radius of pipe wall 

x5 

r_ inner radius of wick 

o 

r„^ hydraulic radius of complete wick structure 
unD 

r^ distance between centers of screen filaments 
f 

resistance of interface at condenser 
iO 
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IE 

\0T 

\e 

^LI 

T 

PCI 

T 

PEI 


8 

AT 

V,. 


resistance o£ liquid vapor interface at evaporator 

resistance of pipe wall at condenser 

resistance of pipe wall at evaporator 

total resistance 

resistance of vapor 

resistance of wick at condenser 

resistance of wick at evaporator 

temperature of outer surface of pipe at condenser 

temperature of outer surface of pipe at evaporator 

temperature of liquid at the liquid vapor interface 
temperature of inner surface of pipe at condenser 

temperature of inner surface of pipe at evaporator 

vapor temperature 
wall thickness of pipe 

thickness of liquid layers between screen layers 



viscosity of vapor 
density of liquid 
density of vapor 

angle within each quadrant over which liquid is present in 
the wick at the evaporator 

thickness of central slab 


The equations developed for computing resistances are given in the 
following list. 

The resistance of circumferential layers of capillary structure and 
working fluid gaps is (see Figure'S) 
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R 

w 


N N-1 

■ E • E 

n=l n=l 


( 1 ) 




R = 
n 


Xg - (n-l)4r^ - (n~l)g 

Xg - nAXj - (n-1) 3 
__ - 

w ax 


( 2 ) 


to 


n,A 


Xg - n4x^ - (n-1) 3 
Xg - n4x^ - n3 

2tt k„ H 
Z ax 


(3) 


Dxy out of the evapoxator suxface is accounted fox by (see Figuxe 4) 


6 = 


2L 

2 


1 + 


"‘MAX 


(4) 


Assume that the active condensex length is il = A , — • 

The wick xesistance at the condensex then becomes 
N N-1 


n=l n=l 


(5) 


R = 
n 


„ x^j - (n-l)4x - (n-1) 3 

W to-^ f-— 

- 4nx^ - (n-1) 3 

2tt k Z ^ 
w cd 


( 6 ) 


R 


nZ 


W ~ ^”^f “ 

Q Xg - 4ni^ - n3 

h \d 


(7) 


The wick xesistance at the evapoxatox is 

N N-1 

I i — 

Q 




1 + 


■^MAX 


2]^n-^X]V. 

n=l n=l 


( 8 ) 
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r - (n-1) 4r - (n-l)e 
- 4n - (n-l)B 

\ = 2pk~JZ 

w E 


(9) 



£n 


- 4nr^ - (n“l)3 

- 4nr^ - nB 

2’ •'j! h 


( 10 ) 


The effective thermal conductivity of the fluid-metal combination in a 
typical single layer is (reference 2) 



The resistance of the pipe wall in the condenser and evaporator sections is 




2n 


■^MAX 


2ir 5. j 
f cd 


( 12 ) 




1 + 


■^MAX 


An 


B 


2ir A„ 
f E 


(13) 


The inter facial resistances at condenser and evaporator ends is approx- 


imated as 
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R 


IC 


Q 4ir r £ , p h. ^ g 
c cd fg ^c 


R 


IE 


1 + 


"^MAX 


( 2.)^/2 ^ 3/2 


2 . 1/2 


4ir r p h- ^ g 
c E *^v fg ®c 


(14) 


(15) 


The resistance of the vapor is computed by assuming fully developed 
laminar flow and accounting for the obstruction presented by the slab 
with a hydraulic diameter. 




8 y A .. 

V eff 


IT p h 
V 




(16) 


’'CHD 2it r^ - 2 5^ + 4 r^ 


(17) 


Hare (reference 1) has developed equations, based on National Bureau 
of Standards data, which give nitrogen properties and stainless steel 
properties in the temperature range of 60“K to 125®K. 

Vapor Pressure [Ib^/ft^] 

= 1.71041 X 10"^ (T)5 - 1.20901 x 10“^ (T)*^ + 3.71275 x 10“^ (T)3 
- 5.70868 X 10(T)^ + 4.28513 x 10^ (T) - 1.25125 x 10^ 

Density of Liquid and Vapor [Ib^/ft^] 

= - 5.8917 X 10~^3 (t) 7 + 4.50297 x lO-^*^ (T)^ - 1.15298 x lO""^ (T)^ 
+ 4.95327 X 10"® (T)*^ + 2.9749 x lO"^ (t) 3 - 5.98552 x 10“l (T)^ 

+ 4.54425 X 10 (T) - 1.21455 x 10^ 


( 19 ) 
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= 1.39324 X 10-13 _ 1.042325 x 10“1° (T)® + 2.638736 x 10"® (T)® 

- 1.14015 X 10“® (T)** - 6.78395 x 10"*^ (T)® + 1.385389 x 10~1 (T)^- 

- 1.07628 X 10 (T) + 3.10045 x 10^ 


Viscosity of Vapor and Liquid [(Ib^ - sec)/ft^] 

= 8.55910 X 10"21 (t)7 - 6.55918 x 10“1® (T)® + 1.70105 x 10"1® (T)® 

- 8.08553 X 10"! (T)*^ - 4.27309 x 10"! 1 (T)® + 8.90377 x 10"® (T)^ 

- 6.k007 X 10"^ (T) + 1.99577 x 10"® 


U- = 4.48282 X 10"! (T)'^ + 2.34251 x 10"! 1 (T)® 

X/ 

- 3.55312 X 10“® (T)2 + 3.14221 x 10"® (T) + 2.16226 x 10"® 


( 22 ) 


Thermal Conductivity of Liquid Nitrogen [ — ] 

k. ■= 1.0970566 x 10"11 (T)® - 9.2427627 x 10“® (T)*^ + 3.090593 x 10"® (T)® 
- 5.1457532 x 10"*^ (T)^ + 4.2210737 x 10“^ (T) - 1.26105 

Heat of Vaporization [Btu/lb ] 


hjg = - 4.11334 X 10"“ (T)® + 2.0908 x 10"® (T)® - 1.43119 x 10"® 
- 1.03235 X 10"® (T)® + 2.61594 x 10"^ (T)^ - 2.40246 x 10 (T) 

+ 8.89614 X 1Q2 


(24) 


Surface Tension [Ib^/ft] 

a = 6.70239 x 10"12 (t)*^ - 4.60497 x 10"® (T)® + 1.19096 x 10"® (T)^ 

- 1.44813 X 10-** (T) + 7.58324 x 10"® 

Ratio of Specific Heats of Nitrogen 

k = 1.572403 X 10"® (T)2 - 8,6844907 x lO"'^ (T) + 1.52913275 (26) 

Thermal Conductivity of Stainless Steel [Btu/ft hr °R] 


k^ = - 4.02016 X 10"® (T)2 + 3.20878 x 10"^ (T) + 1.30266 
where T = temperature in degrees Rankine. 


( 27 ) 
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For any heat pipe of specific structure there are four primary 
variables directly related to the thermal properties of the heat pipe. 
These are T^, T^, T^, and Q. Among these four variables any two are 
independent while the other two are dependent. 

The present numerical procedure involves first specifying Q and 
then making use of the fundamental relationship of temperature, 
resistance, and heat transfer rate to calculate temperatures at several 
locations along the path of heat flow and to calculate the thermal 
resistance of the various sections of the heat pipe. The temperatures 
calculated are only approximate since the simple resistance equation 



is based on the assumption of constant conductivity across the heat 
flow path between the two locations at which the temperatures are 
known. This approximation should not lead to appreciable errors since 
the length of each heat flow path considered is small and consequently 
the value of AT is small. 

The equations which results from applying this method cannot be 
solved analytically because of the presence of high order terms, so 
an iterative solution is necessary. The iteration process used in this 
analysis converges rapidly on the solution, partly as a result of 
knowing the range in which the solution lies, so that the computer 
time required for this portion of the numerical analysis is not large. 
The iterative procedure converged on the solution to an accuracy of 4 
decimal places in an average of only 2 iterations. 

An outline of the numerical procedure is now given 
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I. Specify Q and T 

U 

(Q will vary from 10 Watts to 30 Watts) 

(T^ will range from 60“K to 125°K) 

Approximation: Tjt ^ 

XiX V 

This is an excellent approximation since the resistance of the 
liquid-vapor interface is extremely small. (R^ is typically 
about 0.0001 ®K/Watt) 

The properties of each section will be evaluated at a temperature 
which is the mean of the temperature of the boundaries of the section. 

Algebraic manipulation of the equation Q == 


t3 +( 

4 PCI V 2 2 4 




ct,T, 


ct„T, 


a,T, 


PCI 423 2 


T„a 


C 2 


X 3 2 ^ 2 „ X - ^ 

J-C 2 C ^ 


■3’"C 2tt SL , r^ 

cd B 


- 0 


where ~ ^ 4.02016 x lO"'^ «2 = 3.20878 x 10 


1-2 


) 


PCI 


(28) 


= 1.30266 


II. 

Solve equation 28 for 



III. 

Assume a value for T.^ and compute 

T’ = (T^ + T. 

IV. 

Compute 

from equation 23 

for 

T = T*. 

V. 

Compute kj 

from equation 27 

for 

T = T'. 

VI. 

Compute k 

from equation 11. 




w 

N 



VII. 

Compute 

^ from equation < 

6. 


N-1 

VIII. Compute o equation 7. 

I n* A/ 


Solve Equation 29 for T^ using linear interpolation method, 
repeating steps IV, V, AH) VI. 


IX. 
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X. 

XI. 

XII. 

XIII. 

XIV. 

XV. 


XVI. 

XVII. 


XVIII. 



(29) 


Compute from equation 5. 

Take T' =[Tpgj + 2 . 

Compute from equation 23 at T'. 

Compute k^ from equation 27 at T’. 

Compute from equation 11 at T*. 

Assume values of T and determine T using iteration, repeating 
Jriijl rliX 

steps- XI through XIV and using equation 30, 


T - T 
PEI V 


Q 1 + 


-s-1( 


"‘MAX 


N 


N-1 


-1 = 0 


(30) 






Compute R^- from equation 8. 

Solve equation 31 for T„. 

XIj 

Substituting metal conductivity equation into Q = (T^, - Tpg^)/Rpg gives 


“i / “i “o \ 

— i. T ^ + 1 ^ T + T d — \ T ^ 

4 E \ 4 PEI 2 PEI 2 j E 




2 “2 °1 2 “2 

PEI 2 ^PEI ^ 4 PEI ^ 2 PEI 


Q Jin — 

a T 

p3 _ _2 „ t ^ 

PEI 2 PEI 3 PEI 40 Jl„ 


“3 ) 

)=■ 


E 


(31) 


Compute k „ from equation 27. 
pE 


REPEODTJCIBILrrY OP TEiE 

original page is poor 
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XIX. 

Compute Rpg from equation 13 

XX. 

Compute from equation 24 

XXI. 

Compute p^ from equation 18 j 

XXII. 

Compute from equation 14 

XXIII. 

Compute Rj^ from equation 15 

XXIV. 

^eff = 3-0 ft - ijj/2 

XXV. 

\ 

Compute from equation 21. 

XXVI. 

Compute from equation 18. 

XXVII. 

Compute p from equation 19. 

xxvm. 

Compute from equation 16. 

XXIX. 

Compute RtOT = ^ 


alternate equation 


^OT 


The. complete computer program is given in Appendix A, 

The effects on thermal resistance of changing the size of liquid 
gaps, wall thickness, mesh size, and number of circumferential layers 
of capillary structure has been examined on the computer for a nitrogen 
heat pipe. Table 1 shows some of the geometries studied. Figures 5 
through 21 show how temperature difference and resistance change with 
changes in operating temperature and heat flux for the cases described 
in Table 1/ 


Several figures related to Cases 1 and 6 are included to show 
order of magnitudes for resistances for these two extreme cases. Only 
one plot is given for each of the other cases. It is interesting to 
note that overall temperature difference may increase or decrease with 
increasing vapor temperature at constant heat flux. For example in 
Case 1 the heat pipe temperature difference increases as the vapor 
temperature increases while just the opposite is true in Case 6. 



Case # 

Wall 

Thickness (mm) 

Mesh Size 

Filament 
Radius (mm) 

Distance Between 
Filament Centers 

Number of 
Screen Layers 

Thickness of 
Fluid Layers (mm) 

1 

1.0150 

400 

0.0155 

0.06919 

2 

0.03048 

2 

0.3968 

400 

0.0155 

0.06919 

2 

0.03048 

3 

2.0300 

400 

0.0155 

0.06919 

2 

0.03048 

4 

1.0150 

250 

0.0241 

0.1198 

2 

0.03048 

5 

1.0150 

100 

0.0815 

0.2937 

2 

0.03048 

6 

1.0150 

400 

0.0155 

0.06919 

1 

— 

7 

1.0150 

400 

0.0155 

0.06919 

3 

0.03048 

8 

1.0150 

400 

0.0155 

0.06919 

2 

0.01524 

9 

1.0150 

400 

0.0155 

0.06919 

2 

0.06096 


Table I. 

Values of Parameters for 

Cases Considered in 

This Study 
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The reason for these different trends is that the thermal conductivity 
of liquid nitrogen decreases with increasing temperature while the 
conductivity of stainless steel increases with increase in temperature. 
In Case 1 the conductivity of the metal controls the effective conduc- 
tivity of the fluid metal combination whereas in Case 6, the liquid 
nitrogen controls the effective conductivity of the combination. 



38 


Capillary Limitations 

Writing momentum, energy and continuity equations for steady operation 
of the model heat pipe at capillary limited heat transfer and making the 
standard simplying assumptions the following equation is obtained. 




8y^ il 


IT y. 


Vv 


eff 

TT~ 


(32) 


where 


Q = Capillary limited heat transfer rate 
CL 

CThp Pj 

N = — = "Heat Pipe Number" 

VIt 


a = surface tension of liquid 
hj^ = heat' of vaporization 
p^ = liquid density 
y = liquid dynamic viscosity 

L ^ * 

- pore radius at evaporator surface 

K = — 7 = effective inverse permeability for slab 

based on approach velocity. 

■■ ■ - •+■ ■ — 


n. 




''A 


"B 


K 


A 


= total thickness of slab 
= number of layers of fine mesh in slab 
= number of layers of coarse mesh in slab 
= thickness of a single layer of material A 
= thickness of a single layer of material B 

= inverse permeability for material A based on approach velocity 
= inverse permeability for material B based on approach velocity 
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^eff ~ effective length of liquid path in slab 
b = width of slab 

K = inverse permeability for material at evaporator and 
condenser surfaces based on approach velocity 

L = average distance traveled by liquid in circumferential 
capillary structure at evaporator or condenser (approxi- 
‘ mately 45° arc) 

n^ = number of layers of capillary material on circumference 
6^ = thickness of a single layer of material C 
= axial length of evaporator section 


= axial length condenser section 
= dynamic viscosity of vapor 
= density of vapor 

r^ = hydraulic radius of vapor space 


Approximatfely one hundred different capillary arrangements were 
studied in order, to determine capillary limitations. Table II shows 
geometric parameters for six of the combinations examined. Capillary 
limitations as a function of vapor temperature are shown in Figure 22 
for each of the combinations listed in Table II. 

Transient Analog Computer Studies 

A rather simplified transient model of a cryogenic slab type heat 
pipe with radiator connected is shown in Figure 23. Due to limited 
analog computer capacity relatively- few nodes were used. The equations 
written for this model are; 


Q 


e 


2ir Z 

e 


k- 

_E 


p c V 

(T T) +. PP^P . 
el 2 



In 


(33) 



Table II. Description of Composite Wick Systems Considered in this Study 


Wick 

Composition 

Number 

Screen Mesh Size 

Number of Layers 

Screen Thickness 

- m 

Total Thickness 
of Slab - m 

A 

B 

C 


ng 

’"c 

d.xlO** 

A 

6jXl0 


1 

250 

100 

250 

2 

8 

1 

0.867 

0.314 

0.866 

2.68 

2 

400 

50 

400 

2 

5 

1 

0.744 

0.448 

0.744 

2.39 

3 

400 

30 

400 

2 

4 

1 

0.744 

0.622 

0.744 

2.64 

4 

400 

30 

400 

4 

5 

2 

0.744 

0.622 

0.744 

3.41 


Wick 

Composition 

Number 

Wire Diameter 
"C” Layer 

m X 10*^ 

Pore Radius 
"C** Layer 

m X 10^ 

Inverse 

Permeability-l/m^ 

Effective Permeability - 1/m^ 

K.xlO"^ 

A 

KgXlO"'' 

KgXlO"^ 

n.6. 

K = 6J( ) X 10-7 

1 

0.482 

0.359 

58.4 

2,610 

58.4 

2,716 

2 

0.311 

0.191 

163 

195 

163 

207 

3 

0.311 

0.191 

163 

63.5 

163 

67.3 

4 

0.311 

0.191 

163 

63.5 

163 

69.6 



Q CAP, LIMIT (WATTS) 
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Figure 23 Analog Model 
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Att Z k 
e P 


p pe p w we p A B 


(T, 


Xi) 


47t Jl k 
e w 


tp c V + p c V ] In (r /r ) 
p pe p w we w B C 


(T‘ - T-) 

vap 1 


(34) 


dT 4TT k A 

... yap w_g (T ~ T ) 

d0 [pc V + pc V +2mc]ln (r„/r^) 1 vap 

w we w w wc w a a B C 


4ir k A 

5L_9 CT - T ) 

+ 2m c ] In (r^/r ) '‘^2 ^vap'' 
w we w w wc w a ^ ii o 


(35) 


dT, 


4ir k A 
w c 


(T - T^) 


de [pc V + p c V ] In (r /r_) ^ vap 2^ 

p pc p w wc w B C 


4ir k A^ 

V + p c V ] In (r /r ) ^^c ^ 2 ^ 

p pc p w wc w A B 


(36) 


dT . 4-fr k A 
c _ _E_ c 


4it r A 


(T,, - T ) + 


a c 


de p c V ln(r /r ) ^2 pcVR'r c 

ppcp AB ppcpc 


(T - T ) (37) 


dT 2 tt r A ' e A cr Q 

r _ a c „ s r t ** + 


de me 
r r 


(T - T ) - 


c r m c r me 
r r r r 


(38) 


Limiter 


.(T, - T ) = 
1 vap 


< ^CL 


27t k a 
w e 
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where 


A - area of radiator 
r 

= effective specific heat of adiabatic section and slab 
= specific heat of pipe material 

c^ - specific heat of radiator 

m = mass of adiabatic section and slab 
a 

m = mass of radiator 
r 

= heat flux into evaporator 

Q = heat flux into radiator 

space 

Q = net heat flux from radiator 
r 

R contact resistance between node c and node r 


T = 
e 


vap 

T = 
2 


T - 
c 

T = 
r 

V = 

e p 

V = 
c p 

V 

e w 


temperature node e 
temperature node 1 
- temperature node vap 
temperature node 2 
temperature node c 
temperature node r 
= volume of pipe material in 
= volume of pipe material in 
= volume of wick material in 


V - 




* 1 L ♦ 


evaporator 

condenser 

evaporator 

condenser 


c w 

e = emmissivity of radiator 
- density of pipe material 

= effective density of wick and working fluid 
a = Stefan-Boltzman constant 


6 = time 
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The limiter equation allows one to include the capillary limitation 
in computations. 

Figure 24 shows a schematic of the analog computer circuit. Note 
that inclusion of the radiator introduces non-linear terms in the 
equation. However, these non-linear terms have thus far caused no 
difficulties in the computations. 

Figures 25, 26, and 27 show sample results of some computations 

performed for a nitrogen heat pipe of configuration. Figure 25 shows 

performance for a step change of 5®K in the evaporator temperature. 

Notice that heat transfer at the hot end (Q ) is limited for some time 

e 

due to capillary limitations. The system has essentially stabilized 
after 60 seconds. Figure 26 shows how all parameters vary for a rela- 
tively fast sine wave. Notice that the capillary limitation considerably 
affects heat transfer through the evaporator. As expected, the compu- 
tations indicate progressively smaller oscillations in temperature as 
one moves away from the evaporator and finally the radiator increases 
with time but oscillates very little. There is a considerable phase 
shift between oscillations in different temperatures. Figure 27 shows 
the system changes for a relatively slow variation in evaporator 
temperature. Evaporator heat transfer is again limited by capillary 
restrictions. The amplitude of temperature oscillations tends to be 
more uniform throughout the pipe than in the case where fast oscillations 
were, considered. There are large phase shifts. 

Transient Digital Computer Studies 

The analog computer program described above is limited to small 
transients and thus is of limited use. For this reason a digital computer 
program is now being developed to handle transient computations. The 
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Figure 26. Fast Sine Wave 
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Figure 27 . Slow Sire Wave 
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new approach will allow predictions to be made for various start-up 
situations including start-up from the supercritical state. 

The digital work is now in the early stages of development. 

Fluid d3mamic affects have not yet been incorporated. However, it is 
anticipated that these important affects can be readily included at 
the appropriate stage. 

In the preliminary model now being studied several assumptions 
are made. (See Figure 28) 

Evaporator sadle: lumped mass, contact resistance to pipe wall, 

known heat input or fixed temperature. 

Wall: nodes in r and 6 directions, contact resistance 

to wick. 

Wick; dryout circumferentially as f(Q) innermost layer 

of nodes at same temperature as vapor. 

Vapor: lumped system, includes mass of slab, linear 

temperature drop along tube. 

Adiabatic section 

and condenser: nodes in r,z directions, nodes become active as 

f (Q) , innermost node of wick at same temperature 
as vapor. 

Wick: contact resistance to pipe wall. 

Wall: contact resistance to radiator. 

Radiator: lumped mass with known heat input or fixed temperature. 

Axial conduction: evaporator temperatures averaged for boundary 

condition in adiabatic section, weighted fraction 
of total heat transfer in axial direction subtracted 
from each evaporator node. 




C-n 


Figure 28. Heat Pipe Model 
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Thermal properties: constant at the temperature of the last time step. 

The nomenclature used in the digital approach is; 
area of radiator 


'R 


w 

£ 

c 

A 

e 

“r 

m 


specific heat of radiator 

specific heat of saddle 

specific heat of vapor and slab 

thermal conductivity of xd.ck 

length of adiabatic and condenser section 

length of evaporator 

mass of radiator 

mass of saddle 


s 

m^ mass of vapor and slab 

NIE number radial nodes in evaporator 

NJE' number circumferential nodes in evaporator 
NIC number radial nodes in condenser 
NJC number axial nodes in condenser 
Q heat -input to radiator from space 

heat input to saddle 

contact resistance pipe to saddle or radiator 
r^ inside radius of pipe 

r , radius at node i 


X 

r 

o 


T 

T 


n 

c, , 
n 


E. . 


outside radius of pipe 

temperature of condenser node i,j at time step n 
evaporator node temperature at i,j 


radiator temperature 
K 

T^ saddle temperature 



S3 


vapor temperature 

a thermal diffusivity of pipe wall 

Ax 1/FI ln(r /r^) (see coordinate transformation) 
o 1 

Ay 2ir/NJE (circumferential node width) 

Az A^/NJC (axial node width) 

A 6 time increment 

E emissivity of radiator 

a Stefan-Boitzman constant 

It is convenient to transform from cylindrical to rectangular coordinates. 
In cylindrical coordinates is 

a^T 1 3T 1 a^T a^T 

“aP"^ r ar r ap aP" 

Make the • following substitutions (See Figure 29). 

X ■= In r/r 

c 

y = ^ 

2 = Z 

thus 

ai _ 3T dx , 3T'^ _ 1 3T 

3r 3x dr 3(|> dr r 3x 

a^T 1 

-ap 3r '■ r 3x " 3P r^ 

32 ^ _ a^T 

ap 3y^ 

a^T a^T 

aP = 


1 ^ 

r^ 3x 
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Substitution into V% 

1 1 3T . 1 9T . 1 9T , 9^z 

9x^ ~ 9x 9x 9y^ 9z^ 

1 1 9^1 9^T 

9x^ 9y^ 9z 


Figure 30 shows the transformed computer grid for each of the vari- 
ous regions. The computational procedure for each time step 


A0(n ->■ n+1) is: 


1) Explicit Balance on T 


.. A6 2A0ir r Z 

T + o e 

s “me (NJE)(m )(c )(R 
s s s s c 


j=i 



o e A0 

Vr ~ 


(V) 


Vr 


4) Alternating direction implicit evaluation of evaporator grid, 
(developed next page) ; and 

5) Alternating direction implicit evaluation of condenser grid, 
(similar to 4) 




Becomes 



A 

J.n — an —2- £n — 


Figure 29, ‘ Polar to Rectangular Transformation 
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As an example of the Implicit equations used, consider an interior 
evaporator node for an i sweep. 



The equation can be rewritten as 



References 6 through 9 have been utilized extensively in developing 
the approach described above. 
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CONCLUSIONS 

The primary general goal of this project was to develop techniques 
for predicting transient operation of cryogenic heat pipes. In particular 
the work was aimed towards development of schemes for predicting start up 
from various initial conditions such as those encountered in the super- 
critical regime. In accomplishing these goals it was necessary to first 
study steady state operation. The steady state work included prediction 
of performance limitations and thermal resistances. The transient part 
of the project has teen divided into two areas: subcritical operation 

and supercritical operation. • 

During calendar 1975 the steady state part of the program was essentially 
completed and some results of computations are included in this report. 

The development qf schemes for predicting transient operation is pro- 
gressing well at this time and some preliminary results are included herein. 

It. is significant to note that several graduate and undergraduate 
students have received valuable training while performing tasks under 
this grant. One M.S. thesis, directly related to the project, was pub- 
lished during 1975 and it is expected that another one will be completed 
about July 1976. 
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APPENDIX 


THEKMAL RESISTANCE PROGRAM 




PAGB BLAM.’ NOT fimm 



rjoonnnnnrjn 


PRflGPAM MAIN 


74/70 0PT=1 


FTN 4.4+Roni 


7 


PROGRAM MAIN (INPUT, OUTPUT, TAPF5sINPUT,TAPF6=0UTPUT) 


thermal analysis of a cryogenic heat pipe 

WRITTEN BY OAVTO RUIS 10/29/75 


UNITS OF INPUr QUANTITIES ARE AS FOLLOWS 

TEMPERATURE R Q - — BTLI/HR LENGTH —FT 

OTMENSTON AYTV(?03), AYDT(?03) ,AYQ(203) , AYRT0T(205), AYRPE(?(H), 

C AYR WE (20 5) , A YR IE ( 203 ) , AYRV (20 3 ) , A YRTC ( 203) , A YRWC (203 ) , AYRPC (205) , 
CAYTC(203) 
dimension Ifi(51?) 

DIMENSION TV5(203) ,RT5(203) ,RPE5(?03) , RWE 5 (203) , R TE5 ( 203) ,PV5(203) 
C,RICS(203),RWC5(203) ,RPC5(205) 

DIMENSION rvi (203) , TV1 U203),RT1 (203),RT1 1 (205) ,RWE 1 (203) , 

CRWEl 1 (203) ,Rwri (203) ,RWC1 1 (203) 

RFAL N,KW,KF,KL 

real lca,le,muv,lfff,lcd,lastt 

INTEGER COUNT, CT 
INTEGER sSYM 

Kw(T7V)=Kt(TZV)/((DF/(2.0*RF))*(2,0*(KL (T ZV)/KFf TZV) ) tl)F/(2,0 + Rn 
C-2.ni )+2,n*KL( T7V)/( (DF/(2.0*RF) ) *( (KL( TZV) /KF(TZV) ) *( (2.0*RE }/ 
C(DF-?.0*RF) + j .0) ) )+KL(TZV) / (2.0*RF/(!)F-?.0*RF) + 1 ,0)**2.0 
KF (TYV) =-4.020 16E-5*TYV* *2.0 t3.20B78E-2* lYVfl ,30266 
FUNCH(0TV,I)TAV) = (DTV-TPCI)/( (QMAX/(2,0*PI*KW(DTAV)*LCf)))*SUMA 
C + (i)HAX/(2.0*PT*KL(OTAV)*LCD) )*SUMR1-1 .0 
KL(TXV)=1 .0970S66E-1 t * T X V* * S . 0-9 , 2427627E -9* T X V* *4 , 0 
C+3.090SP3F-6*! XV** 3.0-5. I 457532F-a* T X V * *2 , Of 4 . 22 I 07 3 /E -2* TX V 
C-1. 26105 

FUNCC(OTPFI ,02T AVI = ( (OTPFT-TV)/( (P 1/(2 . "o* THE T A ) ) * ( ( 1 . 0 / ( 2 . 0 *P I 
C*KW(02T AV)*LE) ) *SUMAf ( 1 . 0 / ( 2 , 0* P I *KL ( 02 T A V ) *L E ) ) *SUHH ) ) ) -Q 
FUNCA(Tn=(Al/4.0)*Tl**3.0f(Al*TC/2.0fA2/?.0-Al*Tr/4.0)*Tl** 

C2.0f ( A1 *lC*iC/4.nf A2*TC/2,0f-A3-Al*Tr*TC/?.0-A2*TC/2.0)*T1-Al*Tr. 
r**3,0/4.0-A2*TC* TC/2.0-A3*TC-(QMAX/(2,0*PI*LCD) ) * ALOG (R A/Rbi ) 
FUNCD(inE)=( A1 /«,0)*0TE**5.nt((-Al/4.0)*TPEIf ( A1 /2.0) *TPFTf A2/2.0) 
C*OTF*r)TF+((-A1 /2.0)*TPE I*TPEI-( A2/2. 0) *TPFIf ( A 1 / 4 . 0 ) * TPE 1 * T PE I 
Cf ( A?/2.n)*TPEIf A3)*0TEf ( (-Al/Q.0)*TPEI**3.0-( A2/2.0)*TPEI*TPFT 
C-A3*TPE I-0*ALnr,(RA/Rfi) /(4.0* THETA*LE)) 

CALL PL01S( 10,512,6,00) 
call PL0TMX(34H.fl) 

CALL 21.01(3.0,2.0,-5) 

XP»12.0 
DO 400 L=6,6 
GATE=0.0 

SYMsO 

WRITE (6, 3333) 

REA0(5,71 n)WT,RF,l)F,N,RFTA 
7111 F0RMAT(T7,F12.«, T16,F1 2.8, T3 1 , F 1 2 . 8 , T4 3 , F5 , 1 , T48 , F 1 2 . 8 ) 

WRITE(6, 7222)L 

7222 FORMAT(T70,"CASF «",T77,I2) 

WRITE(6,7777) 

WRTTF(6, 7226) 

7229 FnRMAT(T55,''wT",TS0,"RF",T65,"DF",T7 7,”N",T«7,''BETA'’J 

WRITF.(6, 7230)WT,RF,nF,N,8FTA 


RBPRODUCIBILrry OF -THE 
ORIGINAL PAGF IS POOR 



o n n 


program main 


7A/7'4 


OPT=! 


PTN f»,4»+Ra01 


7?50 F0RMAT(T30,Fl?.ft, T«5, F 15.fl, T60,F1?. 8, T7S, F5. 1 , T82, F 1 2.81 
WRITEfb, 77771 
SL ABTH=R,0F-3 
LEsO.S 
■ urosi.o 
RAsO. 04187 
QS102.39 
iJMAXsl 02.39 
DQs-6.828 
PTS3.141S9 
A1 =-4.02016E-S 
A2s5.20H7flF-2 
A3sl . 30288 
rClNCsO, ^879396985 
RRsRA-WT 

RC =Rrt- 4. ()*RF*N-rt(- TA* (N-1 .0) 

RrHD=fPI*RC*RC-?.0*RC*SLABTH)/(2,0*PT*RC-2.0*SLABTH44.n*RC) 
on 300 K=1,2 

rc=i o«.n 

1)0 190 J=1,200 

1 THFfA=(PI/2.0)/f 1.0+(Q/OMAX)) 

N I = I N r ( N ) 

COMiMfTTNG TPCt AND RPC 

CnUNT = tl 
fcPSl=0.001 
TPCT=TC 
LAST r = TC 

4 FrFltNCAnPCn 
IF(F)S,?0,7 

5 F 1 =F 
TP1=TPCI 

6 TPCl=lPCI+0.5 
F=F0NCAM PCI) 

TF(F)6,?0, 10 

7 F2=F 
IP2=rPCT 

8 TPCT=TPCI+0.S 
FiFlJNCACIPCn 
IF(F)9,20,6 

9 F1=F 
TPlsTPCT 
GO TO 13 

10 F2=F 
TP2=IPCT 
GO TO 1 3 

13 rP3=(TPl*F2.TP?*Ft )/{F2-Fl) 

F5=FUNCA( TP3) 

7A=AHS(LASTT-TP3) 

l.ASTIsTP3 

IF(7A-EPS1 )?1 /21, 1 4 

14 TPtrfP? 

TP2=rP3 
F 1 =F2 
F?=F3 

CnUNT=C0UNT+1 



PROGRAM MAIN 


7 <l/ 7 a OPTsl 


FIN i|.4 + Ril01 


7 


GO TO 

ao NRITE(6, 1030)1PCI 

1030 F0RMAT(T30,”TPCI =", I 38, F7 «P, T50/ "DOUBTFUL" ) 

GO TO 30 
2t IPCI=TP3 

IF{GATE.NF.0,0)GO TO 30 

WRITE (6, 10ilO)TPCI,COUMT,F3 

tOiiO F0RMAT(T30,"IPCT = " , 1 38 , F 7 . 2 , T50 , " COUNT =" , T59, T5, T65, "F =",T69, 
CF8.3) 

WHITE(b,22??) 

30 i Av=nPci+rc)/2.o 

RPC = Al.nG(RA/Rtt)/(2.0*PI*(Al*TAV*TAV4A2*TAV + A3)*LCn*(0/QMAXl ) 

C 

C COMPUriNR TV AND RWC 

C 

. PSUMsn.O 

on /iOT = l ,Ni 
FTsfFLOAT(l) 

TFRM = AL0G( (RB-(F1-1 .U)*0.0*RF-(FI-1 .0)*BETA)/(RB-/4.04FI*RF-(FT-1 ) 
r*RET A> 1 
PSUM=TERM + PSIJM 
4J0 nriNTTMUF 

SllMArPSIH 
PSUM=0.0 
NM1 =N1-1 
on 02T=l,NMt 
FT=FL0AT( I) 

TFRM = ALnG( (RB“/1.0*FI*Pf-(F I-l . 0 ) *BET A ) / { HR- A . 0*F I *RF-F I *BF1 A )) 
PSIif^=TERM+PSliM 
02 CONTINUE 

SUMB = PSIJM 

IFCNHl ,FO.OKSUMR = 0.0 
C 

C ITERATIVE SOLUlfON BY LINEAR INTERPOLATION METHOD 
C 

bO TV=TPCI 

EPS?=o.oni 

bl IAV=(TVtTPCT)/2.0 

S? F=FUNCBn V, T AV) 

IF (F 160,90,70 

60 F1=F 
TVX1=TV' 

61 IV=TV+S.O 
r4V=(TV+TPCI)/2,0 
F=FMNCHI'TV,TAV) 

IF(F)61,90,63 

63 F2sF 

TV2=TV 
GO TO /9 
68 F2=F 

TV2=TV 

70 rv=TV+s.o 

T AV=( TV+TPCT }/2.0 
F=:FUNCB( TV, TAV) 

1F(F)71,90,70 

71 F1=F 
TVXl=TV 



PROGRAM MAIN 


m/TH nPT = t 


fTN «.a + R^IO{ 


7 


7R LnUNIsO 

flO TV3=nvXl-*F2-TV?*Fn/(FP«FtI 

TAV=(TV3+rPCI)/?.0 
cnuNTscnuNT+1 
ZB=ARS(TV3-TV2) 

IF(7B.l F.EPS?)GO TO 91 
F3rFimCB(TV3, TAV) 

F1=F2 
TVX1 =TV? 

F? = F3 
TV2=TV3 
GO TO HO 

90 WRTTEIft, 2010)1 V, COUNT 

2010 FORMAT! rSO/ "TV =" , TS5f Ffi . 2, T65, "COUNT T73, IS, Tft5, "DOUBTFUL " ) 

91 TV=TV3 
IF(GATE.NE.O.O)GO TO OS 
wRnefh,20«(l) IV,C0UNT,F3 

2000 FORMAT(T30,”TV =", T35,FR.?,T50,"COUNT =" , TS8, TS, Tft5, "F =",109, 
CF6.3) 

-^RITF Ct.,2222) 

95 CONTINUE 

RWC = ( 1.0/(2.0*PI*KR(T AV)*tCD*0/0MAX) ) *SU”Al-( 1.0/(?.0*PI*Kl (1AV)*Lf 
C0*0/0MAX> ) *RUMB 
C 

C -AT EVAPORATOR COMPUTATION OF TPE I AND RWF 

C 

TPei=TV 

TAV=(TPEIfIV1/2.n 
F=H)NrCf TPf T, TAV) 
iF(Fnon,i 2 s,ns 
!00 IPlrlPEl 

FI =F 

102 rPF.I = TPFl+S.() 

TAV=( TPEU TV) /2.0 
F = FIINCC ( r Pt T , T AV ) 

IFIFJUjP, 125,105 
lOS TP2=IPEr 

h? = F 

GO TO 1 20 
ns rP2=TPFI 

F? = F 

116 IPET=TPFI+S.O 

F AV=( F-PEI+ rv)/2.(> 

FrFUNt'Cf IPEl , TAV ) 

IF(F) 1 Ih, 125, 1 19 

119 IP1=TPEI 
FI =F 

120 COUNTS 0 
FP83 = 0.00 I 

121 rPi=(TPl *F2-TP2*M )/(F2-F 1 ) 

C0UNT = C0IJNT+1 

TAV=(T.P3+TV)/?,0 

/C=ARS(TP3-TP?) 

TF(7C.LE.EPS3H;0 TO 1?<U 

F3 = FUN(;C( TP3, TAV) 

FI =F2 
IP1=IP2 



OltXGl5>fAT. 



PPObRAM MAIM 


// 4 / 7 /i 


llPTr) 


FTN 


1 


FP-F3 
TP?=TP3 
GO TO l?l 

125 wRITEC6,20B0)7PFTr COUNT 

2080 format ( TSO r '^TPEI , TS7 , F8 . 2 r T67 r **C0UN T =*M75 / 1 S . TH5f "nOUOTFiiL" ). 

120 !F(GATF:.NF,0,0)r,0 TO 130 

126 ^RITE(6,20R0) TP5f COUNT,F 5 

2090 FnRMAT(T50,»TPET =^M 37^ F6 . 2 r T50 , » COUNT = % T58 . T 5 , T 6^4 r " F r",T68, 
CF6.3) 

WRITE(6,2222) 

130 TPET=TP3 

Rwe = (PI/(2.0*THFTA) )*r (1 .0/f2,0*PI*L F ) ) * ( SUM A/K W ( T A V ) +SUMR/KL M AV 

on 

c 

C COMPUTATION OF TF AND RPt 

C 

1 K = TPF I 
F=FUNCD( IF) 

IF(F ) 1 50, IHOr 16S 

150 FF1=TE 

FI rF 

151 TF=TF+l,0 
F=FUNCD( IE) 

TF(F ) 151 / 1B1 , IS^ 

1S3 TF2=:TF 
F2 = F 

GO TO 170 

165 rF2zTF 
F2=:F 

166 TF=TE+1.0 

F = FUNC(M IF) 

IF(FM6M, 181 , 166 
168 TFl = fK 
F1=F 
17U CT-0 

EPS/4 =i0»00 1 

171 TF3=( Tt MF2-TF2^tM )/CF?-F 1) 

CT=Cr+ 1 

/n=ARS(TE3-TF2) 

] F r zc.i F^FHsai on ru iso 

F3 = FUNCi>( TE 

F1=F2 

TEt=TE2 

F2 = F5 

IF2=Tf ^ 

GO TO 171 
IHQ TFrTE3 
F = F3 

rF{OATE.NE*0,0)GO TO 183 
181 WPITK(6, S030) TE.F.cr 

3030 F0RMAT(T50r F8 * 2 , T a5 , , TO 6 , F 1 0 ^ 3 / ThO , » CT =*M6SrlS) 

183 TAV=( TF + TPFn /2.0 

RPE=(PT / l2*0*THFTAn *ALnG(RA/RR)/f 2*0*P'T*LF.*KF(TAV) ) 

C 

C COMPIITAFIUN OF FNTFRFACIAL RESISTANCE 
C 

HFG=/78. l6*(-0* 1 1 33/4E-11 *TV**6.0 + 2.09n8E-8*TV**5»0-l .0 U 1 9F-6* 
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C IV* *4.0-1 .0’^£?35E-^*TV**3.0 + 2.61594E-l*TV*TV-2.a0246E+ t*TV + ft.aP 61 4 
CE + P) 

PV=1 .7I041E-6*1V**S,0-1,20901F-3*TV**4,0+3,71?7SE-1*TV**5.0 
r-S.VOftbfiP +1*TV*TV + 4.?R513E + 3*TV-1 .2512SE + 5 
LCA=LCD*0/RMAX 
C 

C UNITS OF PTC AS COMPUTED HFRF ARF (SFC R / LRF FT) 

C 

RTC=14.40S*TV**2.5/(RC*1.CA*PV*HFG*HFG) 

PIE=(PI/(P.0*THF.TA))*RIC*LCA/LE 

C 

C computation of vapor RESIvSTANCE 
C ONUS ARF (SEC R / LRF FT ) 

C 

LEFF=3.4-LE/2.0-LCD*0/(2.0*QMAX) 

MUV = 8.SSRlE-2t *TV**7.0-6,5591«E-ia*TV**fe.0 + 1 , 70 I OSE- 1 5* I V**5, 0-R . 
C0«b3 3F-14*I V**4.0-4.2 7309E-1 1 * TV** 3 . 0 + 8 . 90 377E-9*T V* I V-f>, 94 007F-7 
C*TV+1 .99S7/F-S 

ROV=l . 39324F-1 3* IV**7, 0-1 , 04P32SE- 1 0*TV**6, 0+P.63B /38F-a* T V** S . 0 
C-1 .1401SE-8*TV**4.0-8. 7839SF.-4*TV**3.0+1 . 30b 3R9E- 1 *T V * T V- 1 .0 76PR 
CE+1*1 V+3. 1 004SF+2 

W0L = -5.«91 7F-1 5*TV**7.0+4.S0297F-10*TV**t).0-l . 1S29ME-7*TV**S. 0 + 4 
rS327E-6*IV**4.0+2,974RE-3*TV**3,0-S.9RbS2F-l*TV*TV+4,S44PSF+l *TV-1 
C.214SSE+5 

RV=8.0*MUV*| KFF*TV* ( 1 .O/ROV-1 ,0/HOL)/(PI*RnV*HFG*HFR*RCHn**4.0) 

C 

C UNIT rOWVFRSION 

C 

SI =0.4097 
kVSI=S1*RV 
RIESI=S1*Hir 
PICS1=S1*RIC • 

S?=l .89R7 
RPCSI=S?*RPr 
PPKSI=S2*RPE 
RWCS J =S?*RWC 
RWES 1=SP*RWF 
S3=0.SS58 
irsi=s3* rc 
rPCTSI=S5*lPtI 
TVSI=S5*IV 
TPEISI=S3*IPF I 
TEST=S3*IF 
051=0*0.2950 

RIOT=RPFSI+RWESI+RIESI+RVSI+RICSI+RWrSI+RPCSI 

S4=304,8 

XRT0T=(TESI-Trsi)/0SI 
DFST=S4*DF 
RFS I=S4*RF 
RFTASI=S4*RFI A 
RASI=S4*RA 
WISI=S4*»JT 
DELI AT=TtSI-irSI 
C 

C WRITING ANP SETTING UP ARRAYS 

C 

IF(K.NE.nGO TO 18b 


BMODUCIBILnY OF THE 
ORIGINAL PAGE IS POOR 



Pi^DGRAM MATN 




opr=i 


FTN /4.<) + R«01 


7 


TVl ( J)=TVSl 
RTlf J)=RTOT 
RWE1 ( J)=RWES1 
RWC1 (J)=RWCSl 
Gn TO 18K 

185 IFCK.NE.I nr,0 TO 187 
TVl 1 ( JlsTVSl 
RT11 (J)=RTOT 
ftWf 1 I (JT=RWESI 
RRCl 1 ( J)=RWCS1 
GO TO 188 

187 IPrK.MF.5)GO TO 186 
?V5(.I)=TVST 
RT5(J)=RTIJT 

RPF5( J)=RPtSI 

RwfcS(.I)=HW£Sr 

R1K5CJ)=HIFS1 

RV5( J)=RVSr 

R1C8(.1)=RTCSI 

RWC5C J)=RWCSI 

RPC5(J)=KPCSl 

188 CONTINUE 
AYTV(J)=I VST 
AYDT (.n=DELTAT 
AYO( J)=QSI 
AYRTOI (.IIsRTOT . 

AYRPKT J)=RPEST 
AYRWFf J)=RNFST 
AYRIEf J)=R IFSl 
AYRV(,n=RVST 
AY«rC(J)=RlCSl 
AYRWC(.n=RwCST 
AYWPrM)=RPCSI 
AYrC(J)=T(,SI 

GA IE = 1 .0 
TC=TC+TCINC 
IQ() CONTINUF 

wRITEf8,33?5) 

3333 FORMAlUHIl 

RPlTE(6,b063) 

8083 F0RMAT(T?0,"I!NITS , T 3S, " TF MPERATURF — K”, T55, "RESIST ANCF -- K 

C/wAT 1 ", T7Q, "HFAT TRANSFER RATE -- WAITS") 

WRITE(b,7777) 

WRIlFffa, 7777) 

WR1TE(6, 7777) 

WRITFI6, 77/7) 

7777 FflRMATdHO) 

WRITE 16,5005) 

500S FORMAT (TS, "Q" , TPR, " T V" , TPO, "OFL T AT" , T38, "R TOT " , T«9 , "RPE" , T 60 , 

C’RWF", T71,’'RIE",T8P,"RV",T93,"RIC",T10i),"RwC",T115,"RPC",T1?8,"TC" 

n 

wRlTEf6,??R?) 

DO POO Jsl ,?00,6 

WRTTE(6,5000 ) AYO(J) , AYTV( J) , A YDT ( J) , A YRTUT ( .1) , AYRPE ( J ) , A Y PWF C ,F) , AY 
CRTEU) , AYRV( J), AYRir.U) , A YRWC ( J ) , A YRPC f J ) , AY T C ( .1 ) 

50 0 0 forma T( TP, F 6.P, T?0,F6.P, T2H,F6.P, T35,E10.0,T/l6,E10.a,T57,E10./J, 
rr66,Et0.a, T79,Kli).4,T90,E10.a, nOl ,F 10.0, T1 12,E10,a, T1P3,F8.P) 
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200 CONTINIJF 
Q=0+D0 
C 

C PLOTTING DT .VS, TV AT CONSTANT 0 

C 

IF(K,NE. 1 )G0 TO 510 

CALL ^SCALE(AYTV,h,0,?00,+l) 

FVT V=AY1 V(201 ) 

0VTV=AYTVC202) 

AYOT( 20 ns 0,0 

CALL SCALFCAYDT, 7.0,201 i+1) 

FVI)T = AYDT ( 202 ) 
nvr>T = AY 0 T( 203 ) 

CALL AXTSC0.O,().O, 17HVAP0R Tfc'MPf HATURF 1 7, 6 .0, 0 , 0 
CALL AXTSIO.O/O.O, 7HDFLTA T, +7, 7 . 0 , RO , 0, FVDT , OVDI) 
510 CUNTINUF 

AY 1 V(?01 IsFVTV 
AYTV{20?)=0VTV 
AYDT(20n=FVf)T 
AYnT(20?)=0V0T 

CALL LTNfcC AYTV, A YO T , 20 0 , 1 , +S0 , SYM ) 

SYM=SYMf 1 
300 CONTINUE 

CALL PLDI f XP, 0.0,-3) 

C ARRAY 35 

c plotting rtot,rwf,rnc, vs. tv 

C 

CALL SCALE! TV 5 , 6.0,200, + n 
FVI V=TVSr?Ol ) 

. 0VTV=TV5(20?) 

RTSC 201 1 = 0.0 

CALL SCALFCRTS, 7.0,201 , 1 ) 

FVR=RT5(202) 
i)VR = RT5(205) 

RT 5 I 201 1 =FVR 
RT 5 ( 2021 =OVR 

CALI AXTS(O,O,0.O, 1 7HVAPUR TEMPFRA T URE , - 1 7 , 6 . 0 , 0 . 0 , 
CALL AXTSfO.0,0.0, IOHTHERMAL RfcS IS T ANCE , + 1 » , 7 . 0 , RO , 
SYN=1 

CALL LINE! TVS, RTS, 200, 1, +50, SYM) 

SYMSSYM+ I 
RWE 5 ( 201 )=FVR 
RwR5(202)=I1VR 
RWC 5 ( 201 )=FVR 
RwC 5 ( 202 )=oVR 

CALL L INEnvS,RWE5,200, 1 , ♦50,SYM) 

5 YM=SYM + 1 

call LTNE(TVS,RWC 5 , 200 , 1 , + 50 , SYM) 

C 

C PLOTTING RPE,RPC, VS. TV 

C 

SYM=SYM+l 

CALL PLOT (XP,0. 0,-3) 

IF(RPES( l)-.GT.RPC5( 1 )) 61 0,615 
610 RPF5(201 )=0.0 

CALL SCALELRPF5, 7.0,201 , 1 ) 

FVR=RPFS(202) 


,FVTV,I)VTV) 


FVTV,nVTV) 

0,FVR,nVRl 
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DVR=HPES(?0^) 

GO TU hl« 

61S RPC^(?On=0.0 

CALL SCALPCRPCS, /.0,P01 » n 
FVR = kPC5(^()?) 

DVR=RPC5(203) 

61H CnNTINUF 

RPES(?01 )=FVR 
RPE'SCPOairOVW 
RPCSCPOl )=FVR 
HPC5f?0?)=f)VR 

CALL AX ISfO.O, 0.0, 17HVAPUP rEMPEKATUPF., - 1 7 , ft . 0 , 0 . 0 , F V T V , 11 V 1 V ) 
CALL AX IS r 0.0, 0.0 , IPhTHFRMAI REST S T ANTE , + 1 « , 7 . 0 , 00 . 0 , F VR , DVR 1 
CALI L IME(1 VS,RPFS,?0n,+1 ,+5rt,SYM) 

SYM = SYM f 1 

CALL LINEd V5,RPC5,?0n, +1 , +SO,SYM) 

r 

C PLOTTING R1E,HTC VS. TV 

C 

SYM=SYMt 1 

CALL PLOI T XP, O.O,- d 
TF (RTF 'i( n ,GI .RTFS (1 ) )650,bS5 
ftSO RTfcS(P0n=0.0 

CALL SLAI FlPlFS, 7.«,P01 , n 
FVR = R7FS< ?0^) 

OVRrRIF'S(POi) 

GO Tf) bS^ 

hSS RICS(?01 )=0.0 

fall scAi F.fPics,7.n,pni , n 

FVP=RICS(?0?) 
i1VR = RirG{?i1 Jl) 
ft'io CONTINUF 

RIEG(?(il )=F VH 
RTFS(POP}=i)VR 
RICSf P01 IrhVR 
Rirs(PO?)=nvR 

call ax ISdI.n, 0.0, 17HVAP0R TFMPF.KATURE , - 1 7 , ft . 0 , 0 . 0 , F V TV , DVTV 1 
CAi.L AXTS(0.(i, 0.0, ISHTriFRMAl RFSISTANCF,4 1 8, / . 0, RO . 0 , F VR, DVP 1 
CALL I INE( 1 VS,RTFft,?00, 1 , +50,SYM) 

SYM = SYM+ 1 

CALI LINE (TVS, RT (.8,800, 1,+ ‘50,3 YM) 

C 

C PLOT UNO RV VS. TV 

C 

SYM = .SYMH 

CALI Pl.ni (XP,0.0.-3T 
CALL SCALKRVS, 7.0,200, f n 
FVR=RVS(201 1 
DVRrRVSf 202) 

CALL AXTSf O.OfO.O, 1 7HVAPf)R T LnPER ATURF , - 1 7 , ft . 0 , 0 . 0 , F V T V , DV T V ) 
CALL AXIS (0.0, 0.0 , JRhTHERmAL REST STANCE , + 1 8 , 7 . 0 , 90 . 0 , FVP , DVR ) 
CALI LIMtnvS,RVS, 200, 1,0,0) 

C ARRAYS a 1 

C 

C rtmt,rwe,rwc 

C 

SYM=SYM+ 1 


ETOODUCIBlLirY OF TEE 
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nPT = l 


PIN 4,a + R4f)l 


CALL PLnTfXP,0.0,-'5) 
call SCALFdVI rf>.0,?00, + l) 

FVTV=TX1 (?01 ) 

DVTV=TV1 (20?) 

KT 1(20 t ) =0.0 

CALL scalecri 1 ,7.o,?oi , n 

FVR=RT 1 (?n?) 

OVR=RT1(?05) 

RT) (?01 )=FVR 
RIU?0?)=nVR 
RWF 1 (P01 )=FVR 
RWK1 (?02)=DVR 
RWCl (?0 1) =FVR 
RWCI (?0?)=OVR 

CALL AX IS (0.0, 0.0, 17H VAPOR TEMPP R A IDHF , - t 7 , (> . 0 , 0 . 0 , F V T V , l)V T V ) 
CALt AXIS (0.0, 0.0, IRHTHEH’^AL RES IS I ANCE , 1 1 8, 7 . 0 , 90 , 0 , FV» , OVR ) 
CALL L INF ( I VI ,RTl,?no, t ,+S0,SYM) 

SYM = SY'>U 1 

CALL L INF f? V 1 ,R'-.EI ,200, 1 ,+50,SYM) 

sym=sym+i 

CA(L LINEflVl, i-'WCl ,200,1, fbO, SYM) 

C 

C ARRAYS tt 1 1 

C 

L WIOl,PWF,RWC 

C 

SYH=SYM+ 1 

CAi I pi.ni (XP,o.n,-3) 

CALI SCALK TV11 0,200, + 1 ) 

FVJ V = 1V11 (201 ) 

i)viv=ivn ( 20 ?) 

RT 1 1 (201 )=0.0 

CAI L SCALP (RI 1 1,7.0,201 , 1 ) 
h VRrRT 1 1 (20?) 

DVRrKIl 1 (?03) 

RT 1 1 (201 )rPVR 
Rin (202)=DVP 
R-wt 1 1 (20n=F VR 
RWEl 1 (202)=0VR 
HWCl 1(201 )=FVP 
RWCt 1 (202)=nVR 

CALI AX IS (0.0, 0.0, IHH thermal RES I S T ANCE , + 1 « , 7 , 0 , 90 . 0 , F VR , DVR ) 
CALL AX IS (0.0, 0.0, 1 7M VAPOR T E MPF R A TURE , - 1 7 , 6 . 0 , 0 . 0 , F V TV , oy T V ) 
CAI I L INF (TVl 1 ,Rni ,200, 1 , +S0,SYM) 

SYM=SYM+1 

CALL LINE n V1 1 ,RWF 1 1 , 200, 1 , +S0, SYM) 

SYM=SYM+ 1 

CALL LINE (TVl 1 ,RwCl 1 , 200, 1 , +SO, SYM) 
call pi m (XP, 0.0,-5) 

400 CONTINUF 

CALI PLOTCl .0, 1 .0,999) 

???? PnRMATdH ) 

99 STOP 

FMD 



